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Abstract. We describe X-ray production in the atmo- 
spheres of hot, early- type stars in the framework of a 
"stochastic shock model". The extended envelope of a 
star is assumed to possess numerous X-ray emitting "hot" 
zones that are produced by shocks and embedded in the 
ambient "cold" medium in dynamical equilibrium. It is 
shown that the apparent lack of X-ray variability on short 
(~ hours) timescales do not contradict a shock model for 
X-ray production. The character of the X-ray variability 
is found to depend on the frequency with which hot zones 
are generated, the cool wind opacity to X-rays, and the 
wind flow parameters, such as mass loss rate and terminal 
speed. 
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1. Introduction 

The X-ray emission from hot stars has proved to be an 
important "window" for investigating mass loss in early 
type stars, and the advent of the latest suite of orbital 
X-ray telescopes are providing new insights. Here we de- 
velop a model to explore the X-ray variability of early 
type stars, first to understand the near absence of sig- 
nificant variability observed so far, and second to deter- 
mine the requisite S/N and time resolutions necessary to 
measure the expected X-ray variability with current and 
future instrumentation. The objects of interest are hot, 
luminous stars with supersonic winds that carry signifi- 
cant mass, such as OB and Wolf-Rayet (WR) stars. The 
main parameters characterising these stellar winds are the 
terminal velocity (Uoo) and mass loss rate (M). Typical O 
star mass loss rates are M « 10~ 7 Af Q yr~ 1 with terminal 
speeds woo ~ 1000 - 3000 km s _1 . For Wolf-Rayet (WR) 
stars M is larger at around 10 _4 ' 5 Af Q yr _1 but with ter- 
minal speeds like those of O stars. Such strong outflow re- 
sults in a predominance of anomalously strong and broad 
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emission lines in the spectra of WR stars. The high mass 
loss rate during the WR phase has a significant effect on 
its evolution and the chemical enrichment and mechanical 
energy input to the interstellar medium. 

It is widely accepted now that stellar winds of OB stars 
are maintained by radiation pressure in numerous lines of 
heavy elements (Castor et al. 1975; Pauldrach et al. 1986). 
But the acceleration of high mass loss WR winds is though 
to be due to multi-line scattering of photons (Lucy & Ab- 
bott 1993; Springmann 1994; Gayley et al. 1995). Non- 
LTE, moving atmosphere modelling codes have been de- 
veloped (e.g., Hamann & Koesterke 2000; Hillier & Miller 
1998), which adopt radiative and statistical equilibrium 
and assume a monotonic velocity law a priori. In the 
framework of the so-called "standard model" the emis- 
sion lines are formed in a spherically symmetric, time- 
independent, dense, smooth stellar wind that is photoion- 
ized by the hot core. However, there has long been growing 
evidence that at least some assumptions of the standard 
model are oversimplified. A pointed example is the X-ray 
emission from hot star winds, first detected as discrete X- 
ray sources with Einstein (0.2-4.0 keV) (Harnden et al. 
1979; Seward et al. 1979; Seward & Chlebowski 1982). 

Early UV observations of hot stars revealed the ap- 
pearance of superionisation, especially the presence of Ovi 
which would not be expected in the winds of hot stars. 
Cassinelli & Olson (1979) suggested that Auger ionisation 
by X-rays in the context of a coronal model could explain 
the existence of these highly ionised species. However, the 
hypothesis of a hot corona in the outer parts of a stellar 
atmosphere was ruled out due to the fact that the soft X- 
ray flux is not sufficiently absorbed at low energies (e.g., 
Cassinelli & Swank 1983). According to a phenomenolog- 
ical model proposed by Lucy & White (1980) and devel- 
oped by Lucy (1982), dynamical instabilities should arise 
in line driven stellar winds leading to shock generation 
throughout the flow. Cassinelli & Swank (1983) found that 
the shock model of Lucy (1982) produced too low an X- 
ray flux and also that the X-rays were too soft to ex- 
plain the Einstein observations of the O stars. They also 
pointed out that the spherically symmetric shock model 
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would lead to strong variability of the X-rays, which was 
not seen. To explain the lack of clear evidence for vari- 
ability of the X-rays from hot stars they proposed that 
the X-rays are not from spherically symmetric shocks but 
instead from shock fragments in the wind. Rosat (0.2-2.4 
keV) observations have brought new insight in the prop- 
erties of X-ray emission from hot stars. Kudritzki et al. 
(1996) have found interesting relations between the X-ray 
emission and the wind properties from an analysis of 42 
O star spectra. For example, they confirm the relation 
between X-ray luminosity and stellar bolometric luminos- 
ity as Lx ~ 10~ 7 £bo1 that was previously known, and 
also find that Lx scales with the wind kinetic luminosity 
0.5Mti^. For the WR stars, the data is of lower quality 
and consists almost entirely of passband fluxes. However, 
based on Einstein observations, Pollock (1987) was able 
to find that on average, N-rich WR stars (WN) tend to 
be more luminous than the C-rich WR stars (WC), pos- 
sibly attributable to the difference in abundances of the 
two types. The Rosat data have also revealed that, un- 
like their progenitors the O stars, the X-ray luminosities 
of single WR stars are not correlated with LboIi wind mo- 
mentum Mvoo, wind kinetic luminosity 0.5Mti^, or WR 
subtype (Wessolowski 1996; Ignace & Oskinova 1999). Re- 
cent data from Chandra and XMM-Newton of OB 
stars (Schulz et al. 2000; Kahn et al. 2001; Waldron & 
Cassinelli 2001) and WR stars (Maeda & Tsuboi 1999) 
are just becoming available. These newer data, mostly of 
O stars, are showing rich X-ray line spectra, the analysis 
of which are yielding some unexpected results. We refer 
the reader to the above mentioned papers for a descrip- 
tion of these data and here concentrate on the topic of 
X-ray variability in single stars, a topic that has not much 
been addressed from a theoretical perspective. 

Dynamical instabilities in line-driven winds have been 
studied extensively for OB stars (Owocki & Rybicki 1984; 
Owocki et al. 1988; Rybicki et al. 1990). Model computa- 
tions predict shock velocity-jumps ranging from 500 km 
s _1 to 700 km s _1 , implying post-shock temperatures and 
emission measures that are marginally adequate to explain 
the observed thermal X-ray emission (e.g. Hillicr et al. 
1993, Feldmeier et al. 1997a). Gayley & Owocki (1995) 
have shown that similar to OB stars, the instability mech- 
anism should also operate in the denser WR winds. 

In terms of observed X-ray variability, Crowther & 
Willis (1996) have presented a review of available data 
based on Rosat observations. For ( Ori, Berghofcr 
& Schmitt (1994) claimed to detect X-ray variability 
over a 2-day period, reaching about 30% in hard X- 
rays (0.64-2.38 keV), while the soft band data (0.16- 
0.5includcgraphicsl keV) remained about constant. For ( 
Pup, Berghofer et al. (1996) found evidence for a modu- 
lation in the harder (0.9-2.0 keV) X-ray emission and in 
the optical data with a 16.7 hour period, and suggested 
that this reflects periodic changes in the base wind den- 
sity of the 04 f star. On the other hand, the study of 



a Ori showed no significant X-ray variability. Except for 
WR6 (HD 50896) and WR 1 (HD 4004), no X-ray variabil- 
ity studies for single WR stars have been undertaken. For 
WR 6, Willis & Stevens (1996) reported variability at the 
< 30% level on timescale of < 1 day, together with larger 
epoch changes. No significant spectral shape changes were 
found. In the case of WR 1, no evidence of significant X- 
ray variability was reported by Wessolowski & Niedziclski 
(1996). Thus far, no short-time scale X-ray variability on 
the order of the wind flow time (i?*/voo ~ hours) has been 
observed in O and WR stars. 

In this paper we consider the application of stochastic 
inhomogeneous wind models, as widely accepted necessary 
components to explain spectral variations at optical and 
UV wavelengths, to the X-ray emission from single O and 
WR type stars. We restrict our analysis to shock models of 
X-ray formation in hot stars and do not consider alterna- 
tive models (see discussion in Waldron & Cassinelli 2001). 
In Sect. 2, we describe a stellar wind model consisting of 
numerous zones filled by hot X-ray emitting gas. Two fol- 
lowing sections use two different approximations of the 
optical depths of the cool material, which absorbs X-ray 
emission. In Sect. 3 we apply the so-called exospheric ap- 
proximation treating the X-ray emitting zones as spherical 
shells. In Sect. 4 we generalise our model by using formal 
integration of the radiative transfer equation for randomly 
distributed wind shocks. We investigate the crucial pa- 
rameters of the variability and their relation to the input 
parameters of the model. A summary discussion of the 
results is presented in Sect. 5. 

2. Stochastic model for the distribution of X-ray 
emitting material. 

We consider a spherically symmetric and time- 
independent stellar wind that is a mix of "cool" and "hot" 
gas in dynamical equilibrium. The minor hot gas compo- 
nent gives rise to X-ray emission. We assume that hot gas 
is present in the form of spatially separated compact vol- 
umes characterised by temperature Tx and density which 
arc different from the temperature and density of the am- 
bient cool wind gas. 

It is common in the literature to attribute wind 
"clumps" to dense formations with temperatures much 
lower then the temperatures of the X-ray emitting mate- 
rial, sometimes referred to as DWEEs (discrete wind emis- 
sion elements, Lepine & Moffat 1999). These clumps are 
presumably responsible for observed line profile variabil- 
ity in the optical spectra of WR stars, or the formation of 
DACs (discrete absorption components) observed in the 
optical and UV spectra of some O stars (Brown et al. 
1995). Such cold clumps and hot X-ray emitting zones may 
well have related origins. Feldmeier et al. (1997a) consid- 
ered the formation of X-ray emitting zones as a conse- 
quence of shell collisions due to the propagation of reverse 
shocks in the stellar wind. In their hydrodynamical sim- 
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ulations they demonstrated the existence of dense shells 
moving according to the stationary wind velocity law and 
small-mass, high-velocity blobs. Copious X-rays are pro- 
duced when a fast blob rams into an outer shell. These 
X-rays thus originate over a fairly narrow range in radius. 
In this respect, recent observations (Schulz et al. 2000) 
seem to confirm such a scenario, except that the expected 
strong variability is not observed. 

The general description of a source with a continuous 
spatial distribution of temperature (e.g., a radiative cool- 
ing zone behind the front of a shock wave) is given by the 
differential emission measure which includes integration 
over all surfaces with the same temperature, such as over 
multiple independent shocks in a stellar wind (Craig & 
Brown 1976). But it is necessary to note that an overall 
consistent picture of dynamical, structured, thick stellar 
winds is still far from being complete. Although detailed 
models do not predict the X-ray sources to be isothermal 
with unique temperature Tx amongst them, the isother- 
mal approximation for the overall wind structure should 
be fairly close to the actual case because of the short radia- 
tive cooling times in these dense stellar winds (Feldmeier 
et al. 1997b). The two-temperature approximation with 
isothermal X-ray emitting zones embedded in a cooler out- 
flowing wind is able to reasonably well reproduce the ob- 
served X-ray spectra from hot stars (Hillier et al. 1993) 
and might be used to place constraints on physical pa- 
rameters such as the volume filling factor (e.g., Oskinova 
et al. 2001). 

Therefore, we will consider a spherically symmetric 
and, on average, time independent outflow that is a homo- 
geneous mix of "cool" and "hot" gas in dynamical equi- 
librium. For isothermal optically thin hot X-ray zones, 
the particular choice of Tx = 10 7 K tends to maximise 
the X-ray emission in the ROS AT ba,nd (see Ignace et al. 
2000). We use the Raymond-Smith cooling function (Ars) 
for optically thin plasmas (Raymond & Smith 1977). For 
WR stars, a rough correction of Ars for non-solar abun- 
dances is given in Ignace et al. (2000). The ambient cool 
medium (< 10 5 K) is as described by the standard model 
and is optically thick to X-rays. 

In the two-temperature approximation, the ratio of 
emission measures of hot and cool components of the wind, 
EMx/EM w , is proportional to the ratio of volumes filled 
by hot (Vx) and cold (V w ) material and to the square of 
their density ratio D s . Thus, the "filling factor" is 



EMj 



(i) 



where /x™ and /j£ are the mean mass per free electron of 
the cool or standard wind component. The X-ray emission 
arises primarily from two-body processes-collisional exci- 
tation and radiative recombination. Therefore, the volume 



emissivity j v (energy per unit time per unit mass) is pro- 
portional to the square of the density, 



= fxpl 



A„(T X ) 
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where fif~ is the mean molecular weight per ion of the hot 
gas and mn is the proton mass. A„ (ergcm 3 s _1 ) is the 
cooling function at the energies E concerned and p w is 
the density of the standard wind is 



Pw(r) 
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By the definition of filling factor in Eq. ([!]) , the total spe- 
cific luminosity (erg s _1 ) emerging from the wind 



Lx(E) = / x A 1/ (Tx)EM ¥ 



(4) 



This then is the monochromatic X-ray luminosity at any 
instant in time, since our hypothesis is essentially that 
the emission measure of X-ray emitting material (or cor- 
respondingly the filling factor) is time variable. 

It was pointed out by Cassinelli et al. (1996) that the 
X-rays from a wind shock will be seen only after that shock 
has moved into a region where the cool wind attenuation 
is optically thin to the X-rays. As we shall show, the wind 
is quite opaque at most X-ray energies up until the flow 
has reached a substantial fraction (;> 50%) of its termi- 
nal speed. Thus, we can assume that the X-rays emerge 
primarily from distances where v(r) ~ v^. For the wind 
attenuation, K-shell absorption by metals in the cool wind 
is the dominant opacity source. Assuming photo-electric 
absorption by K-shell electrons 



k w (E) = d v (E)l 'pYrriH 



(5) 



where pf is the mean molecular weight per ion of the 
cool wind and oj (E) is the photo-absorption cross-section. 
(Actually, the appropriate molecular weight to use is that 
per nucleus, but in hot star winds, there is no neutral gas, 
and so every nucleus is an ion.) The ratio n-Jni defines the 
relative abundance by number for atomic species j, where 
rij is the number density of ions. The X-ray absorption 
(apart from prominent K-shell edges) can be represented 
roughly by a power-law in energy: k w (E) « kqE~ 7 , with 
7 in the range 2-3, depending on ionisation structure and 
chemical composition (see, e.g. Hillier et al. 1993; Cohen 
et al. 1996). The parameter kq is a constant that depends 
sensitively on the abundances. 

To proceed, we now develop a phenomenological model 
in which the wind inhomogcncitics develop at random 
times near some inner boundary Rq. Then they propa- 
gate radially according to a monotonically increasing ve- 
locity law v(r). Let us introduce dimensionless notations 
for velocity u — v(r)/voo, and for distance x = r/Rg. 
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Further, we normalise time to the characteristic time 
7f = Ro I Vqo , which is equal to the flow time in the specific 
case i?o = R*- Then for unitless time we have t = T/Tf. 

To avoid confusion we would like to clarify here, that 
further out, we shall use two descriptions for the distri- 
bution of the hot gas. In the exospheric approximation 
described in Sect. 3, we assume that the hot material is 
present in the form of spherical shells. The exospheric ap- 
proach with spherical shells are simplifying assumptions, 
but the model does allow us to gain some insight into the 
nature of X-ray variability from wind shocks. In Sect. 4, 
we allow for randomly distributed wind shocks of arbi- 
trary geometrical form (as long as the hot zones are not 
exceedingly "large", to be defined later), and we consider 
the full radiative transfer problem for the emergence of 
X-rays from the wind flow. 

For the motion of the hot zones (e.g., shells), we con- 
sider an ensemble of zones filled by hot gas propagating 
in the radial direction with constant velocity u = 1. The 
different zones in our model are labelled by their times of 
appearance t at the distance xq = 1. The first zone is at 
distance xq — 1 at the moment of time to — 0. At time 
tfc when shell number k is crossing xq = 1, the first zone 
is at distance x 1 (t^) = xq + tfc with tfc = Yli=i ^» an d 
zone number j is at distance x-^tfc) = xq + (tfc — tj). Here, 
we choose to use the time intervals <5tj as drawn from uni- 
formly random numbers in the range [0.5, 1.5] multiplied 
by a dimensionless parameter (T): 



St; 



"P([0.5, 1.5]) ■ (T), 



(6) 



and so (tj) = (T). The parameter (T) is in fact an average 
separation in time between the appearance of subsequent 
zones and is measured in flow times 7f. Further we will 
explore the influence of the average time separation (T) 
on the character of the variability. 

The parameter (T) is introduced because detailed 
models of X-ray production demand the presence of an ini- 
tial perturbation of the stellar wind. For example, in mod- 
els which depict X-ray production mechanisms in O stars 
(Feldmeier et al. 1997b), the existence of large scale tur- 
bulence near the wind base is postulated. As it will be- 
come clear below, the average time separation (T) is the 
only essential parameter of the models under considera- 
tion. We thus wish to emphasis that in principle, one can 
expect (T) to be an observable value, which could be in- 
ferred from the X-ray emission light curve analysis as the 
average separation between the peaks of intensity. 

The emission measure of each zone is <SEM = n^nf^SV. 
Using Eq. (|3|), the squared number density inversely scales 



with the fourth power of distance: n c n[ 
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where 



(7) 



On the other hand, the volume of a zone SV with fixed 
radial thickness and solid angle grows with the square of 



distance SV oc x 2 . This implies that the contribution of a 
zone to the total X-ray luminosity determined by its emis- 
sion measure SEM has an overall decline that goes as x~ 2 . 
This is a plausible treatment when mass conservation for 
the hot shocked material is assumed so that the continuity 
equation Eq. (||) can be applied. 

3. Exospheric model for the time dependent 
X-ray emission of single stars 

3.1. Modelling time dependent X-ray emission 

Let us assume that the observed X-ray emission arising 
from the hot gas emerges only from radii exterior to the 
optical depth unity surface of radius n (exosphere), with 
X-rays at smaller radii completely attenuated. In this ap- 
proximation the X-ray luminosity has the proper scaling 
with mass loss and other key wind parameters. The ex- 
ospheric radius r\ of optical depth unity is (Owocki & 
Cohen 1999) 



n {E) 



M 
47rw 0< 



(8) 



The emission measure of the cool material outside the ex- 
osphere which is assumed to be transparent to the X-rays 
is given by 



EM, 



4-7T 



n> e (l-W(n))r 2 dr, 



(9) 



where W(t*i) is a dilution factor defined to be w ri /47r 
where u) ri is the solid angle subtended by the exosphere 
and the parenthetical term accounts for geometric occul- 
tation by the exospheric surface of radius r% . 

Let us make the simplifying assumption that X-ray 
emitting material is present in the form of spherical shells. 
All shells have the same small fractional thickness C = 
Ai? s hoii/-Ro at a given radius in the envelope. Recall that 
we assume isothermal X-ray emitting zones and thus the 
parameter C is not to be misconstrued as a cooling length 
but rather the geometrical fractional thicknessof a shell 
filled with the hot gas. 

In our phenomenological exospheric model, X-ray 
emitting shells are seen only beyond radius r% — x\Rq, 
so that effectively xq = x\. The shells start their motion 
outwards in the atmosphere at quasi-random times tfc. At 
that moment, the X-ray flux (6 Fx) is increased by the new 
source. Then owing to the decrease in emission measure 
of the zone with distance (<5EMx oc x~ 2 for constant shell 
thickness), the incremental flux SFx decreases. Thus, the 
total X-ray output of an ensemble of X-ray emitting zones 
will be observed to vary with time. 

The radius x\ has a strong dependence on mass-loss 
rate, chemical composition and energy (see Fig. 1). This 
means that in the model under consideration, the X-ray 
emitting zones will cross the surface of optical depth unity 
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-0.2 0.0 0.2 
Log(E) [keV] 

Fig. 1. Radius of optical depth unity as a function of X- 
ray energy. Dashed-dotted line is for typical parameters 
of 01 stars, the solid line is for typical WN stars and the 
dashed line indicates WC stars, which are enhanced in C 
and O but deficient in N. Note the strong K shell edges 
due to CNO elements 



at different distances for different energies. Therefore their 
emission measures at the moment when they appear for 
an external observer will be different. Obviously then, one 
may expect the character of variability to change as a func- 
tion of wavelength and for stars of different abundances. 

The emission measure of an individual shell j at time 
t fc is 



<5EM,- 



47ri?^<< • C [xi(E) + (t fe - tj)Y 



(10) 



Combining Eqs. (0) and (|l0|), the emission measure of the 
ensemble of hot shells becomes includegraphicsinclude- 
graphics 



EMx(£,t fema J 



4irRlnl£ ■ 



fcmax 

E 

k=l 



(11) 



Here the integral over the volume is replaced by summa- 
tion over all X-ray emitting shells outside the sphere of 
x\(E), and fc max is the maximum number of shells under 
consideration in the wind. The value of fc max is chosen so 
that the most distant shells make a negligibly small con- 
tribution to the emission measure (i.e., 8EMx{x m a,yi) — > 
for x > ii). We choose fcmax ~ 10 5 which allows us to 
account for the emission of a shell till it reaches a distance 

■^max fcmax (-0 ^ 10 ■ 

The X-ray luminosity with time is Lx{E,t) = 
EMx(E, t)A„(Tx). Given that hot gas zones are a func- 
tion of radius only, and not of latitude or azimuth, a rela- 
tively simple expression for the relative variability, a^/Lx, 
can be derived. Recognising that the X-ray emission from 
any given zone drops rather steeply with radius, both 



the emission and variability will be predominantly gov- 
erned by zones that appear at the x\ radius. Hence the 
fractional variability can be estimated based on Poisson 
statistics. The number density of zones above a given ra- 
dius x scales as l/ccf; the volume in that vicinity scales 
as a;?; so the number of zones near x% is N oc x%. Conse- 



quently the relative variability scales as a^/Lx oc 1/y/xi. 
Given that x\ oc k v {E), we can expect the relative vari- 
ability to be a function of energy. For purposes of illustra- 
tion, if k u (E) oc _E~ 7 , then the relative variability will be 
&l/Lx oc E^l 2 , which is an increasing function of energy. 
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Fig. 2. Variations of emission measure with time for O, 
WN and WC stars. The solid line represents the vari- 
ation near 1.5 keV foincludegraphicsr a star with an 
O star chemical composition. The dash-dotted line is for 
a WN star. The dotted line is for a WC star. Note that 
(T)=10 for these models. From the numerical simulations 
the ratios of standard deviation a to the average value 
of emission measure are <r/(EM) = 1.0 for the O star, 
ct/(EM) = 0.38 for the WN star and a /(EM) = 0.13 for 
the WC star. 



Figure 2 shows changes in the emission measure of the 
shell ensemble, computed using Eq. ([ll]). We took the val- 
ues of x\ corresponding to 1.5 keV for O, WN and WC 
stars. As is clear from the figure, the character of vari- 
ability depends on the chemical composition. As can be 
seen from Fig.^, the radius of optical depth equal unity 
near 1 keV for a WC star is located about an order of 
magnitude further away in the wind than is the case for 
a WN star. This leads to an overall reduction of the X- 
ray emission, and some suppression of the variability in 
WC stars at this energy. In Sect. 4 we revisit these simu- 
lations using the formal solution of radiative transfer with 
angle dependent optical depth, instead of the simplified 
exospheric approximation. But first we make some com- 
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Table 1. Typical values of C for WR 114 (WC5) and WR 1 
(WN5) Q 





WR114: 


WR1: 


(T) 


fx = 0.175 


fx = 0.0831 


0.1 


6.2 ■ 10~ 3 


5.1 ■ 10" 3 


1.0 


4.5 • 10~ 2 


2.5 ■ 10~ 2 


10.0 


4.4 • 10" 1 


2.0 ■ 10" 1 



a Values of fx and stellar parameters are from Ignace 
et al. (2000), Hamann & Koesterke (1998), and Koesterke & 
Hamann (1995). 



ments on model constraints obtainable from filling factor 
data. 

3. 2. Using filling factors to constrain model parameters 

In Ignace et al. (2000), lower limits to the filling factors 
of most putatively single Galactic WR stars were deter- 
mined from an analysis of Rosat All-Sky Survey (RASS) 
observations. A natural question arises: Is it possible to 
extract information about the spatial distribution of hot, 
X-ray emitting gas by knowing fx from RASS passband 
observations in the framework of the model under consid- 
eration? Let us again assume that the fractional thickness 
of the X-rayemitting layers C does not depend on distance. 
By definition of Eq. (0) and using Eq. (|p) and Eq. ([Tl|) , we 
have: 



fx = CD 2 S 



E 1 



l dE 



(12) 



Here, we neglect the occultation term in Eq. ( J9Q. Note 
that X\ — xi(E) is implicit. One infers from Eq. (|l^) that 
in fact there are two dependent parameters in the prob- 
lem: (T) and C. Let's assume that it is possible to infer 
(T) from observations of the variability of X-ray emission. 
Then making assumptions about the density ratio D s , one 
may determine £ if Lx (equivalently fx) is known. There- 
fore from Eq. , we may estimate C from known values 
of fx- In Tablejlf we give typical values of C using D s = 4 
(the same value as in Baum et al. 1992 and Hillier et al. 
1993) for WN and WC stars. 

So far, we have used the exospheric approximation to 
demonstrate that the X-ray variability is strongly influ- 
enced by the opaque cool material of the stellar wind and 
therefore (a) is a function of energy and (b) depends on 
the chemical composition of the stellar wind (see Fig.^J). 
Also, it was shown that lower limits to filling factors can 
be used to infer lower limits to the thickness of spherical 
X-ray emitting shells (as illustrated in Tab.[|), if the aver- 
age time interval between subsequent shocks (T) is known 
from observations. 



4. Formal solution for X-ray emitting zones 

While it appears that the majority of WR and O stars 
display basically spherical winds on the large scale, there 
is compelling observational evidence that the winds are 
clumped on small scales (e.g. Lepine & Moffat 1999). In 
Lepine & Moffat (1999) a phenomenological model was 
introduced, wherein WR winds are depicted as consisting 
of a large number of randomly distributed, radially prop- 
agating, discrete wind emission elements. This model was 
used to simulate line profile variability patterns in emis- 
sion lines from a clumped wind. We adopt this approxi- 
mation for the X-ray producing region and consider the 
X-ray emission as arising from analogous optically thin 
zones with temperature Tx, which as before are embed- 
ded in the cooler bulk wind flow. However, the zones are 
no longer treated as spherical shells. 

We treat these hot zones as piecewise elements, pre- 
sumably produced by shocks. They expand adiabatically 
moving radially according to a /3-law for the velocity. 
Zones start their motion at randomly distributed moments 
of time from random but uniformly distributed positions 
at an imaginary spherical surface of radius Rq - the min- 
imal predicted radius for shocks to occur. 

The main difference with the exospheric model from 
the previous section is that the optical depth of cool mate- 
rial now depends not only on radius but also on the spher- 
ical polar angle as seen by an observer. Further we now 
compute a proper line-of-sight optical depth to each hot 
emitting zone, in contrast to the exospheric approximation 
for which hot zones are either completely attenuated, or 
completely unattenuated. To simplify the numerical sim- 
ulation, we do require that the emission measure of every 
zone <5EM is the same at a given radius and that the spa- 
tial size of each zone is small enough to assume the wind 
attenuation to every point in the zone is approximately 
constant (e.g., the zones cannot be hemispherical shells, 
since the wind attenuation would strongly vary across such 
a structure). 

In this section we adopt the usual cylindrical coordi- 
nate system (e.g. Lamers et al. 1987). The coordinates 
are: distance from the star x, impact parameter for the 
line-of-sight p, and distance along the line of sight z, all 
normalised to the stellar radius. The angle between the 
line of sight and the radial vector is 9 = arccos /i, with 
fj, = 1 if the radial vector points to the observer. For any 
point in the wind, x, p and z are related by 



z = fix — n(p 2 + z 2 ) 1 / 2 



(13) 



With our assumptions, the emission measure changes 
with distance from the star as 5EM X = 5EM (l/x) 2 , with 
5EMo the emission measure of a zone at distance xq . The 
luminosity of each zone for which the optical depth of the 
ambient wind is r w is 



SL V = A,(Tx)6EM x e- 



(14) 
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Table 2. Representative optical depths to an X-ray emit- 
ting zone propagating radially with 9 = 45° assuming a 
typical WN star a chemical composition. 



U [v(r)/ Deo] 


0.1 


0.9 


1.0 


E [keV] 




0.5 


3.4 ■ 10 3 


3.5 ■ 10 2 


3.5 ■ 10 1 


1.0 


5.5 ■ 10 2 


5.4- 10 1 


5.5 


1.5 


2.1 • 10 2 


2.0 ■ 10 1 


1.9 



a With Voa = 2500 km s" 1 and M = 10 _4 - 4 M 



©yr" 



and the luminosity of a hot gas zone occulted by the stellar 
core is assumed to equal zero. Applying the equation of 
continuity Eq. (|) with Eq. (|), we can derive the wind 
optical depth to be 



a(E)M 



r dz_ [°° dz 

I 2 ( \ = T 0\ E ) / 2 1 N >( 15 ) 

J z x z u(x) J z x z u(x) 



where x = \/ p 2 + z 2 . Now the optical depth for each zone 
is angle dependent. For soft (0.4-2.4 keV) X-ray energies 
at radii where the wind acceleration takes place, the opti- 
cal depth is generally quite large (see Tab.||). 

Although our numerical calculations can account for 
the X-ray emission from zones at small radii, Table |^ in- 
dicates that their contribution to the total X-ray flux will 
be quite negligible owing to the rather extreme attenua- 
tion. So, we simplify our calculations by considering only 
those zones which have already reached the constant ve- 
locity region in their motion throughout the wind, hence 
u = 1. Using Eq. (|l5|), it is possible to obtain an analytic 
expression for the optical depth (MacFarlane et al. 1991): 



ME) 



dz 

2 



to(E) 



sin ( 



(16) 



So, the optical depth depends on the location of a given 
zone in the wind as well as on the energy. 

In Fig. |^ the change of X-ray flux for a single X-ray 
emitting zone along with its motion through the flow is 
shown. The luminosity grows rapidly with distance owing 
to the decrease in optical depth by the overlying wind. In 
general the zones which are propagating in directions close 
to cost? = 1 have higher luminosity due to the fact that 
the optical depth is smaller in this direction. The emission 
rises to a peak followed by a decline. At this point, the 
attenuation is negligible, and the decreasing emission is a 
consequence of the inverse square fall-off of the emission 
measure with distance. The distance x pea k at which an X- 
ray emitting zone yields peak emission depends strongly 
on the angle 9. Further, noting that c5EMx ~ x~ 2 and 
t w ~ a; -1 it can be shown (from an analysis of Eq. (|i~4|)) 
that the maximum of the curve occurs when t w (E, 9) = 2. 

Let us assume that a stellar wind contains numerous 
X-ray emitting zones. Clearly an ensemble of such zones 




Fig. 3. Change in monochromatic X-ray luminosity Lx at 
IkeV (solid line) and 1.2 keV (dashed line) of a single X- 
ray emitting zone with distance X in a WN star wind. The 
zone moves along a trajectory at 9 = 45° in the envelope. 

will also produce time variable X-ray emission. The rise 
in luminosity of a single zone is much sharper than its de- 
crease after peaking. Suppose that at a given time and a 
given energy, the X-ray luminosity of only one zone could 
reach its maximum. Subsequently, the emission from this 
zone will be decreasing, and until the next zone reaches its 
maximum, the emission of the whole ensemble will drop. 
Thus, in the regime of constant expansion the separation 
between maxima of emission reflects the average time sep- 
aration (T) . We may expect larger levels of variability on 
a short time scale (order of flow times) for larger values 
of C. So, we may conclude that the observed lack of vari- 
ability suggests values of (T) ~ 1 and therefore values 
of C ~ 10~ 2 , which are marginally consistent with esti- 
mations of cooling lengths (e.g. Hillier et al. 1993). It is 
necessary to point out here that the available observations 
so far have not been capable of detecting such small fluc- 
tuations of X-ray flux. 

To proceed further, let us recall that the total X-ray 
luminosity depends on the random variable tj. For the an- 
gular distribution of hot zones, we select cos 9 as a uniform 
random variable in interval [—1, 1], and for the azimuth, 
4>k is uniform random in the interval [0, 2tt]. It should be 
noted that although cosflj. is uniform, r Wj fe is not. The to- 
tal emission for an ensemble of zones is a summation over 
all contributors with the appropriate attenuation: 



T tot 



(E) = A V (T X )SEM ]T 



1 



(1 + t k ) 2 



(17) 



The wind will be more opaque to X-rays at soft energies 
than to those at hard energies. Therefore, we expect the 
phcnomcnological picture wherein an X-ray emitting zone 
has peak flux at higher energy earlier than it has peak flux 
at softer energy. 
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Fig. 4. Change in narrow band X-ray luminosity Lx at 
1 kcV with time for WR1 (WN5) (upper curve) and for 
£ Pup (04If) (lower curve). The average separation be- 
tween X-ray emitting zones (T) grows from the left panel 
to the right (T) = 0.1, 1, 10. The flow time for the 
WN type star is if « 0.15 h and for the O star is if « 0.9 h. 

Using £EMo as a parameter of the model, we avoid 
direct references to the density of hot material or volume 
of the zones filled by this material. In our model at each 
moment of time, only one X-ray emitting zone is at a 
given distance from the inner boundary Rq. This is similar 
to the propagation of subsequent spherical shells which 
was described in the previous section in the framework 
of the exospheric approximation. In Section 3.2 we have 
shown that the fractional thickness of a spherical layer C 
and (T) are not independent parameters but coupled with 
each other (see Eq.|lJ). To place constraints on <5EMo, let 
us assume now that SEMo is the same as the emission 
measure of a spherical shell with thickness C located at 
distance Rq with corresponding density no as in Eq. (Q) . 

The number of hot gas zones and their initial emission 
measure <5EMo is set by the average separation in time 
(T) . The assumption of the constant filling factor leads to 
the correlation between average time interval and initial 
emission measures of hot zones. This means that small hot 
zones with small separation in time should produce small 
changes in X-ray flux. On the other hand, large sized hot 
zones produce changes in X-ray flux with a more signif- 
icant amplitude. To complete the picture, it is necessary 
to bear in mind that the degree of variability also depends 
on chemical composition. 

Fig. ^ presents numerical simulations of short-scale 
time variability of the monochromatic X-ray luminosity 
at IkeV for two stars with different cool wind opacities, 
namely O stars and WN stars. Clearly, the character of 
variability strongly depends on the average time separa- 
tion (T) and the spectral type of the star. The mass loss 
rates of WN stars are much higher then those of O stars; 
therefore, we attribute the difference in attenuation be- 



tween WR stars and O stars as the main cause for differ- 
ent character of variability shown in Fig. ^ owing to the 
fact that x\ is much bigger for WN stars than for O stars. 

The changes in luminosity seen in Fig.[| are, in fact, 
due to the superposition of light curves for many zones 
propagating in different directions similar to those shown 
in Fig. Q The average separation between subpeaks on the 
curves from Fig.^ reflects the average time between two 
subsequent zones having peaks in their luminosity. How- 
ever, when (T) is small, the flux fluctuations are negligi- 
ble (order ~ few per cent) and simply cannot be resolved. 
The interesting question to address is what is the plau- 
sible range for the parameter (T). Observed fluctuations 
in X-ray emission for £ Ori (04f) (increase in the count 
rate of « 30% for 2 days) and ( Pup (09.51a) (modula- 
tions with a period of 16.7 hours and amplitude < 10%, 
Berghofer & Schmitt 1994, Berghofer et al. 1996) suggest 
rather small values of (T) (< 1 flow time), assuming that 
the X-rays form in wind shocks. 

As shown by our numerical simulations, on the time 
scale of several hours, the variability of X-ray luminosity 
might be negligible for stars of spectral type WN. The 
reasons for this apparent lack of variability are that in 
addition to small size and short time separation of hot 
zones, the total optical depth for the indicated energies is 
quite large in the WN star winds. As can be seen from 
Eq. (|r^), the exponential term e~ T ™^ suppresses the difer- 
ence in values of Lx ■ Although not shown in Fig. ||, our 
modelling reveals that the X-ray luminosity of WC stars 
is almost constant due to their quite opaque stellar winds. 
As obvious from the figure, even the transparent winds of 
O stars may not demonstrate detectable levels of variabil- 
ity in X-rays on short time scales of about 1 hour in the 
case of small (T) . 

Although on a time-scale of hours, the X-ray variabil- 
ity may be quite small, the dynamical flow time r/v(r) 
for distances r > 100i?* is some 10 5 seconds, of order 
a day. So we have performed simulations of changes in 
flux over longer time intervals. The results of this calcu- 
lation are shown in Fig.^j. Surprisingly, even for (T) = 1, 
the variability becomes significant (~ 10%) on time scales 
of thousands of hours. In this case, the separation be- 
tween peaks in the light curves do not reflect the time 
separation between X-ray emitting zones. To understand 
the character of this variability let us consider Eq. (|l7|) . 
The time dependent luminosity is governed by two uni- 
formly distributed random variables Sti , and cos 9i . Let 
us assume, for simplicity, that zones are launched within 
constant time intervals equal to unity and concentrate on 
random variations in r w . Then 8ti = 1 and from Eq. ([Tt]) 

Elmo*'-* ' (18) 

k 

From Eq. (|l6|), T Wl k oc To#fc//csin#fc and cos 0k is a random 
variable from [—1, 1]. Thus, when the number of realisa- 
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^ 4 



1.5 keV [Fluxx 100] 



1.0 keV 



1.5 keV [Fluxx 150 (WC)] 
1.0 keV (O) 



12 3 

Time (in 1000 flow times) 

Fig. 5. Change in monochromatic X-ray flux at IkeV, 
1.5 keV and 3keV with time for WR1 (WN5) (solid lines), 
for ( Pup (051a) (lowest curve) and a typical WC star 
(dashed line). The average separation between X-ray emit- 
ting zones is (T) = 1. 



tions k is large enough, even for rather opaque winds with 
To 3> 1, the variance of the flux is significant. 

The likelihood of detecting variability of X-ray flux 
increases drastically at hard energies. This is due to the 
strong dependence of the optical depth on energy. To il- 
lustrate this, if k w (E) ~ E~~< then using Eq. ( |l8| ) and 
neglecting the second order terms, the standard deviation 
of the X-ray luminosity becomes ctl ~ E 2 (M/woo)~ 2 . 
This is an increasing function of energy and decreasing 
function of the stellar wind density. As numerical simu- 
lations show, in the case of random times t;, the general 
trend of the dependence of variance on energy and density 
is an increasing function of energy. The only condition is 
that the variance of emission measures of the hot zones 
should be smaller than the variance of the optical depths 
associated with them. 

Fig.^ represents the ratio of the standard deviation 
ctl to the average value of the X-ray luminosity as a func- 
tion of energy, computed using Eq. (17). Apart from ion- 
isation edges, Fig.^J clearly shows that the relative vari- 
ability is an increasing function of energy. Fig. || confirms 
the basic scaling of variability with parameter x\ that 
was obtained under the exospheric approximation with 
&l/Lx oc 1/y/xi. Presumably, this will fail for large val- 
ues of (T) ~ x\ , because then there are relatively few hot 
zones in the outflow. At very low (T) one hits another 
regime where the flow approaches homogeneity. 

As seen in Fig. p] for (T) = 1, the amplitude of variabil- 
ity for O stars can reach as much as 80% near 1.5 keV but 
only 10% at softer energies around 0.5 keV. At the same 
time, the relative variability is only a few percent near 1 
keV for the WN stars. For WC stars, the level of X-ray 



variability is negligible. Therefore, it is not surprising that 
with such low sensitivity X-ray emission detectors such 
as Rosat PSPC (0.2-2.4 keV) or Einstein IPC (0.2-4.0 
keV) , the detected X-ray flux of WR stars appeared to be 
fairly constant X-ray sources. 




0.0 0.1 

Log (E) [keVl 

Fig. 6. Dependence of the ratio of standard deviation ctl 
to the value of the X-ray luminosity (Lx) versus energy for 
typical abundances of O stars (solid lines) and for the av- 
erage chemical composition of a typical WN star (dashed 
line). 



5. Summary and conclusions 

We present simulations of the expected X-ray variability 
for early-type stars in the framework of a shock model 
for X-ray production. We assumed that the optically thin 
hot X-ray emitting material is embedded in a cool X-ray 
absorbing stellar wind, described by the standard model. 
For such a medium we used the concept of filling factor, 
which is the ratio of emission measures of X-ray emitting 
and absorbing material. Assuming an isothermal X-ray 
emitting hot gas, the filling factor is proportional to the 
ratio of volumes filled by hot and cool components and to 
the square of the ratio of densities of hot and cool mate- 
rial. To examine the basic processes of X-ray emission and 
absorption we considered two different models. 

To begin, we employed the exospheric approximation 
with angle independent optical depth. In this approxima- 
tion we used a model of the stellar wind where hot gas 
is present in the form of spherical shells propagating with 
constant velocity. The emission of such an outflow is time 
dependent. This allowed us to derive a expression restrict- 
ing the fractional thickness of a spherical shell emitting 
X-rays by using filling factors. The lower limits on filling 
factors for isothermal X-ray emitting material were de- 
rived in Ignace et al. (2000). Using these data, we found 
lower limits to the thickness of the spherical shells filled by 
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hot gas to produce the observed level of X-ray emission. 
It was shown that the thickness of an isothermal spherical 
shell can be inferred from the analysis of the X-ray light 
curve and the observed filling factor. 

Further, the model of an envelope consisting of a 
number of individual radially propagating X-ray emitting 
zones was developed by analogy with stochastic wind mod- 
els. The total X-ray luminosity of such an atmosphere 
was obtained by using the formal solution of the radia- 
tion transfer equation for angle dependent optical depth. 
We did not place any constraints on the form, volume or 
density of the zones of hot gas except to demand that the 
total emission measure of each zone is the same as the 
emission measure of a spherical shell at the same distance 
from the stellar core, and that the zone must be rela- 
tively "small". This allows us to restrict the number of 
independent parameters when performing the numerical 
simulations. The rationale for doing this was that broad- 
band X-ray luminosities are known for most hot stars. 
Therefore, any combination of parameters describing the 
distribution of hot gas should provide the observed level of 
X-ray emission for each particular star as its first priority. 

The two major simplifying assumptions of our mod- 
elling are those of isothermal X-ray sources and smooth 
cool material. We may speculate that including radiatively 
cooling sources will lead to a decrease in the level of vari- 
ability. On the other hand, taking into account the clumpy 
structure of the background wind will lead to an increase 
of the variability of X-ray output. A detailed consideration 
of these effects are a matter of future work. 

Despite its appoximative nature, the analysis de- 
scribed here revealed the basic properties of variability 
of X-ray emission for early type stars. Some of the main 
points are: 

1. The apparent lack of short time-scale variability (order 
of an hour) of X-ray emission cannot be considered as 
a deficiency of shock models for X-ray production. 

2. We may expect stochastic variability on long time- 
scales (thousands of flow times), especially for optically 
thick winds. 

3. The level of X-ray variability depends on chemical 
composition and density of the stellar wind and differs 
for stars of different spectral classes even for similar 
mechanisms of X-ray production. It is governed by the 
opacity of the cool material and is substantially lower 
for the more opaque winds of WR stars and particu- 
larly WC stars. 

4. The dependence of the wind opacity with energy means 
that the X-ray emission may be highly variable in the 
part of the spectrum, where the cool material is op- 
tically thin and practically constant in optically thick 
parts. Therefore the hard energies are specially apt for 
detection of variability of X-ray emission. Recall that 
existing claims of X-ray variability show differences be- 



tween soft and hard passbands (e.g., possibly detected 
for ( Ori by Berghofer & Schmitt 1994). 
5. Whether X-ray variability is detected or not for a given 
energy provides valuable information about the spatial 
distribution and properties of X-ray emitting material, 
and if variability is detected, its dependence on energy 
would be especially telling of the wind structure. 

Acknowledgements. This research was supported by a PPARC 
grant (LMO, JCB) and NASA grants NAG5-9964 (RI) and 
NAG5-9226 (JPC). We acknowledge discussions with K. Gay- 
ley, M. Hendry and S. Rauzy, and we thank the referee 
L. Koesterke for comments beneficial to this paper. 

References 

Baum E., Hamann W.-R., Koesterke L., et al. 1992, A&A 266, 
402 

Berghofer T.W., Schmitt J.H.M.M., 1994, ApSS 221, 309 
Berghofer T.W., Baade D., Schmitt J.H.M.M., et al. 1996, 
A&A 306, 899 

Brown J.C. Richardson L.L., Antokhin I., et al. 1995, A&A 
295, 725 

Cassinelli J.R, Olson G.L., 1979, ApJ 229, 403 
Cassinelli J.R, Swank J.H., 1983, ApJ 271, 681 
Cassinelli JR., Cohen D.H., Sanders W.T., 1996, In: Tanan- 
baum H., White N., Sullivan P. (eds.), Proc. Throughput 
X-ray Spectroscopy workshop, Harvard Press, 294 
Castor J.I., Abbott D.C., Klein R.I., 1975, ApJ 195, 157 
Gayley K.G., Owocki SR., Cranmer SR., 1995, ApJ 442, 296 
Gayley K.G., Owocki SR., 1995, ApJ 446, 801 
Crowther P., Willis A.J. 1996, In: Vreux J.-M., Detal A., 
Fraipont-Caro D., et al. (eds.), Proc. 33 r d Liege Int. As- 
troph. Coll., Wolf-Rayet Stars in the Framework of Stellar 
Evolution. Liege, Univ. de Liege, 109 
Cohen D.H., Cooper R.G., Macfarlane J.J., et al. 1996 ApJ 
460, 506 

Craig I.J.D., Brown J.C, 1976, A&A 49, 239 
Hamann W.-R., Koesterke L., 1998, A&A 335, 1003 
Hamann W.-R., Koesterke L., 2000, A&A 360, 647 
Harnden F.R., Jr., Branduardi G., Gorenstein P., et al. 1979, 
ApJ 234, L51 

Hillier D.J., Kudritzki R.P., Pauldrach A.W., et al. 1993, A&A 
276, 117 

Hillier D.J., Miller D.L., 1998, ApJ 496, 407 
Ignace R, Oskinova L.M., 1999, A&A 348, 45L 
Ignace R, Oskinova L.M., Foullon C, 2000, MNRAS 318, 214 
Feldmeier A., Kudrizki R.P., Palsa R., et al. 1997a, A&A 320, 
899 

Feldmeier A., Puis J., Pauldrach A.W.A, 1997b, A&A 322, 878 
Kahn S.M., Leutenegger M.A., Cottam J. et al. 2001, A&A 
365L, 312 

Koesterke L., Hamann W.-R., 1995, A&A 299, 503 
Kudritzki R.P., Palsa R., Feldmeier A., et al. 1996, In: Zimmer- 
man H.U., Triimper J., Yorke H. (eds.), Rontgenstrahlung 
from the Universe, MPE Report 263, 9 
Lamers H.J.G.L.M., Morton D.C.,1976, ApJS 32, 715L 
Lamers H. J. G. L. M., Cermti-Sola M., Perinotto M., 1987, 

ApJ 314, 726 
Lepine S., Moffat A.F.J., 1999, ApJ 514, 909 
Long K.S., White R.L., 1980 ApJ 239, 65L 



Oskinova et al.: X-ray variability modelling 



11 



Lucy L.B., White R.L., 1980, ApJ 241, 300 
Lucy L.B., 1982, ApJ 255, 286 
Lucy L.B., Abbott D.C., 1993, ApJ 412, 771 
MacFarlane J.J., Cassinelli J. P., Welsh B.Y., et al. 1991, ApJ 
380, 564 

Maeda Y., Tsuboi Y., 1999, BAAS 31, 1541 

Oskinova L.M., Brown J.C., Cassinelli J. P., et al. 2001, In: 
Moffat A.F.J., St. -Louis N. (eds.), Interacting Winds from 
Massive Stars. ASP Conference Series (in press) 

Owocki S.P., Castor J.I., Rybicki G.B., 1988, ApJ 335, 914 

Owocki S.P., Cohen D.H., 1999, ApJ 520, 833 

Owocki S.P., Rybicki G.B., 1984 ApJ 284, 337 

Pauldrach A.W.A, Puis J., Kudritzki R.P. 1986, A&A 164, 86 

Raymond J.C., Smith B.W., 1977, ApJS 35, 419 

Pollock A.M.T., 1987, ApJ 320, 28 

Rybicki G.B., Owocki S.P., Castor J.I., 1990, ApJ 349, 274 
Schulz N.S., Canizares OR., Huenemoerder D., et al. 2000, 
ApJ 545, L135 

Seward F.D., Forman, W.R. Giaccononi, R. Griffiths, et al. 

1979, ApJ 234, L55 
Seward F.D., Chlebowski T., 1982, ApJ 256, 530 
Springmann, U., 1994, A&A 289, 505 
Waldron W.L., Cassinelli J.P., 2001, ApJL 548L, 45 
Wessolowski U., Niedzielski A., 1996, MPE Report 263, 73 
Wessolowski U., 1996, MPE Report 263, 75 
Willis A.J., Stevens I.R., 1996, A&A 310, 577 



